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Abstract 

We show that, for a fairly large class of reversible, one-dimensional cellular au- 
tomata, the set of additive invariants exhibits an algebraic structure. More precisely, 
if / and g are one-dimensional, reversible cellular automata of the kind considered by 
Takesue PQ, we show that there is a binary operation on these automata V such that 
VK/) ^ ^(/ V g), where ip(f) denotes the set of additive invariants of / and C denotes 
the inclusion relation between real subspaces. 

1 Conserved quantities and symmetries 

Additive invariants for cellular automata (CA) have been widely studied over the last 
two decades, especially in connection with the fundamental role they play in physical 
modeling. The problem of finding explicit conditions for a one-dimensional CA to have 
conserved densities of a given size (in a sense that will be described in this article) has 
been solved by Hattori and Takesue in their remarkable 1991 paper [2]. Interestingly, 
they expressed those conditions in terms of discrete current laws, thereby recovering a 
fundamental property of continuous systems in a discrete context. 

This striking analogy with continuous models suggested that a discrete version 
of Noether's theorem - an elegant and fundamental connection between conserved 
quantities and symmetries of physical systems - might be formulable. On the other 
hand, the proof of Noether's theorem is based on the properties of differential operators 
which are absent in a discrete context. 

For these reasons we decided to approach this problem first from a computational 
point of view. We considered a class Tj of second-order, reversible CA based on 
a binary alphabet defined over a fixed "neighbor window" / and, for each / £ fj, 
computed all its additive conserved quantities, the set of which may be identified with 
a finite-dimensional vector space denoted by ip(f)- Then we looked at the equivalence 
classes of CA rules given by the relation 
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f~g ip(f) = ip(g) 

and observed that, if we choose two representatives / and g from the same equivalence 
class ip(f), then there is a pointwise binary operation V on Ti such that 

W)^(/V 9 ). (1) 

This relation is interesting because it relates a dynamical system to its invariants. Since 
Noether's theorem endeavors to relate the latter with the symmetries of the system, 
we thought it was worth investigating this problem further. In this work we present a 
proof of Eq.([T]) for a large class of reversible, second-order dynamical systems defined 
on an arbitrary finite set. 

2 Physical motivations 

In a seminal paper [I], Takesue argues that a family of fully-discrete models called 
Elementary Reversible Cellular Automata (ERCA) can be used to study the thermo- 
dynamic behavior of large dynamical systems. ERCA are a family of one-dimensional 
reversible CA with two Boolean independent variables at each site. It is possible to 
associate to these systems certain additive conserved quantities which can be regarded 
as a form of energy. Also, discreteness of the variables implies that the phace space 
volume is preserved under the dynamics, and therefore the statistical mechanics of the 
model can be constructed. 

Examples of physical applications of these quantities are given, for instance, by 
Pomeau [6] and Boccara [31 3] who have used them successfully to model traffic. In 
most of these works, special symmetries of the rules are considered and conjectures are 
drawn on how they relate to the associated set of additive invariants. 

In our approach, we give primacy to the equivalence relation which identifies two 
rules if they have the same set of invariants. The main result is a theorem which shows 
that symmetries of the rules are reflected in an algebraic structure on their respective 
equivalence classes. Since "energies" are hydrodynamic quantities which one hopes 
to see emerging from the dynamics, our result can be seen as a contribution toward 
understanding how certain qualities of the dynamics reflect on the symmetries of the 
associated physically relevant quantities. In the next section we start the discussion 
by generalizing ERCA to an arbitrary, finite alphabet. 

3 Generalized second-order, reversible CA 

The instantaneous state of many discrete dynamical systems is most conveniently de- 
scribed by n dependent variables, each of which takes its values in a set S, with 
\S\ = N < oo. We denote the j'th dependent variable at time step t by Xj(t) £ S, for 
< j < n and t > 0, so that the state of the system at time t 

x(t) = {x (t),xi(t) } . . .,X n -l(t)} 

may be thought of as an n- vector. CA dynamical systems of this kind endowed with a 
shift-invariant map [5]. 
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More precisely, consider the set Ip = {—(3, —[3 + 1, . . . , 0, . . . , (3 — 1, /?} ((3 £ Z) and 
a function 

f:S T — ► S. 

Then a one-dimensional cellular automaton is the dynamical system defined by the 
triple 

(S,P,f), 

with the dynamics being described by the following map: 

Xi(t + 1) = f(Xi-p(t), ... , Xi(t), x i+ p(t)). 
By defining xf (t) := {xi_p(t), . . . , Xi(t), . . . , Xi + p(t)}, the above can be rewritten as 

Xi (t + l) = f(x?{t)). 

We refer to the associate global map F : x(t) i— * x(t + 1) as the global map induced by 
/• 

In what follows, we will also assume periodic boundary conditions, which means 
that X(i+j) = X( i+ j) m0f } n and we refer to the system so defined as a cellular automaton 
(CA) of first order and neighborhood size (3. 

We now define the concept of second-order reversible CA, mentioned in the Intro- 
duction. Let 

7 : S 2 — ► S 

be a second-order reversible map; that is 7 satisfies 

Vx,y€s ■ -y(y,i(y,x)) = x. 

Then we can define 

f:S z f>xS — > S 

fa(t + l),Xi(t))»7(f(xi(t + i)),Xi(t))=Xi(t + 2). (2) 

Consistent with the notation introduced above, we call / a second-order CA of neigh- 
borhood size (3 and reversible combiner 7. Since the binary version of these objects 
was first introduced in [2], we will also refer to them as Takesue-type cellular automata. 
When it will be clear from the context, we will drop the tilde from the notation. 
We notice that: 

f(x i (t + l),Xi(t + 2))=7(f(x i (t + i)),x i (t + 2)) = 

= j(f(x t (t + 1)), 7 (/(xi(i + 1)), Xi(t))) = Xi (t). (3) 

Hence, by Eqs.© and ([3]), we deduce that the second-order dynamical system F in- 
duced by /, 

F(x(t + 1), x(t)) = x(t + 2) 

implies 

F(x{t + l),x(t + 2)) =x(t) 
demonstrating that F is indeed a reversible dynamical system. 
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4 Additive conserved quantities 



Consider a second-order reversible CA / of neighborhood size (5 on a periodic lattice 
of size N. We define an "energy density function" 

e : S 2a+1 xS — > R 

(xf(t) ! x i (t-l))^e(^(t),x i (i-l)). 

Remark 1 In what follows, two distinct neighborhoods will appear: One is relative 
to the automaton (neighborhood size (3), and one is relative to the conserved density 
(neighborhood size a). However, the a-neighborhood will appear exclusively as an argu- 
ment of e, whereas the f3 -neighborhood will appear only as a property of the automaton. 
Hence, in order to simplify the notation, we can drop the superscripts a and (5 with no 
ambiguity. 

We define the total energy En of the system at time t associated to e as: 

JV-l 

E N (x(t),x(t - 1), e) = <Xi(t), Xi (t - 1)). (4) 

i=0 

When the density e is given, or clear from the context, it will be dropped from the 
arguments for E. 

We say that e is a additive conserved density for f if, for any t > and iV > we 
have 

E N {x(t),x(t-l),e) = E N (F{x(t),x{t-l)),x(t),e). (5) 

According to this definition, in order for e to be a conserved quantity, an infinite number 
of conditions must hold: namely, conserved quantities are solutions to the system of 
equations defined by Eq. ([5]) for any N G N + . In their seminal paper [2] Hattori and 
Takesue proved that only a finite number of linear conditions are sufficient to find all 
the solutions to Eq. ([5]). 

These conditions can be rephrased in a slightly different, though equivalent, fashion. 
Indeed, let us denote by i^(f,N,a) the set of all the additive conserved densities of 
neighborhood size a for / on a periodic lattice of size N. Then we have: 

Lemma 1 If there is Nq such that ip(f, Nq, a) = ip(f, Nq + 1, a) then 

4>(f,M,a) = (f,N ,a), VM > N 

Lemma 2 For any f and any a there exists Nq < oo such that 

^(f,N ,a)=^(f,N + l,a). 

For the proof of these lemmas we refer to j2]. Let us notice that the above result 
tells us that the set of invariants for a given automaton / is, for N sufficiently large, 
independent of the size of the periodic lattice on which the system evolves. Therefore 
we can drop the dependence on N and define the set of additive conserved quantities 
for / to be ip(f,a) = ip(f,N ,a). 
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5 An algebraic property of additive invariants 



In this section we prove that additive invariants for the second-order, reversible CA 
considered above possess a natural equivalence class structure. 

Remark 2 The choice for the space-time window neighborhood (in particular, the 
choice of the "central" point of the neighborhood at time t — 2) which defines the den- 
sity e is crucial to our analysis. The idea is to force the window for e to have the 
same "shape" as the one for the second-order automaton f : In the next section the 
consequences of this choice on our argument will appear more clearly. 

Using the notation introduced above, let us consider a second-order reversible CA 
with S = {0, 1}, (3 = 1 and a = 1. This corresponds to the following time evolution: 

Xi (t + 2) = f(xi(t + l))®Xi(t) 

where Xi(t) G {0,1} for all i, s, and © is addition modulo 2. Since we set a = 1, 
the density e has the form e = e(xj_i(t), Xi(t), Xi+i(t), Xi(t — 1)) and can therefore be 
represented by a vector with 16 components. 

Let us now fix the window: That is, let us denote by J- the set of all the functions 

{0,1} 3 ^{0,1}, 

and, for any / £ J, let </>(/) be the set of all second-order additive conserved densities 
for / with a = 1. The equivalence relation 
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then induces a partition of T into equivalence classes. We did an exhaustive computer 
search and found that ~ partitions J- into 21 distinct equivalence classes. 

Example 1 The following independent row vectors, written in matrix form as: 



( 1 

1 

-1 

1 


V o o 





1 




1 



2 1 

D 1 

-2 

I 

-1 -1 



■1-10 



2 00 -1 00001\ 

1 1 1 

-10 1 10 

00 -1 01000 

-10 1 10 

11 00000/ 



are all the possible conserved quantities of the CA represented, using Wolfram's nota- 
tion |7|/, by 8,64 an d 72. Therefore {8,64,72} comprises a class in J 7 / ~. 

In the above example, CA 8, 64 and 72 correspond, respectively, to the following 
binary sequences: 

8^(0,0,0,0,1,0,0,0), 64^(0,0,0,1,0,0,0,0), 72^(0,0,0,1,1,0,0,0) 

If we define / V g to be the bitwise inclusive or of the binary representations of / 
and g, we can easily verify that the class formed by 8, 64, 72 is closed with respect to 
this operation. Interestingly, this closure property holds for all the classes except two. 
However, the following slightly weaker property holds in this context: 
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Conjecture 1 If f,g belongs to the same equivalence class, then we have that 



W)^(/V ff ). 

This is interesting, since it says that there is an "invariant-preserving" operation defined 
on the class of CA with a given window function. In the next section we will present 
a lemma which will be used to prove, in Section [71 a generalization of this result to an 
arbitrary finite alphabet. 

6 The main lemma 

Following the notation introduced in the previous section, let / be a second order 
reversible cellular automaton defined by 

/ : S 2/3+1 xS — > S 

(Xi(t), Xi (t - 1)) ^ j{f(xi(t)), Xi {t - 1)), 

where 

7 : S 2 — ► S 

is such that 

Vx,y€s ■ -y(y,i(y,x)) = x. 

Notice that, on a binary alphabet, the only non-trivial function with this property is 
addition modulo 2 (up to conjugation), which we denote by ©. Notice also that, since 
© is associative and commutative, it satisfies 

(x (y z)) = (y (x © z)). 

The next lemma shows that these two properties of the reversible combiner are sufficient 
to prove Conjecture [3 

Lemma 3 Suppose that f is a second-order reversible cellular automaton and 

2. ry(x,j(y,z)) = j(y,j{x,z)). 
Then, if e is an additive conserved quantity density for f and g we have that 

e[si(t),7(/(si(i)),Xi(t-l))] =e[x i {t), 1 {g{x i {t)),x i {t-l))] (6) 
for any i 6 {0, 1, . . . , N - 1}, t > 1. 

Proof. Since the system is reversible it is enough to prove the claim for case t = 1. 

First notice that, if /(£j(l)) = g(xi(l)), the claim is trivial. Then, let x(0) = 
{xo(0), xi(0), . . . , xtv_i(0)} and x(l) = {xo(l), . . . ,xjv-i(1)} be the initial con- 
ditions. 

We set Xf(2) := F(x(l),x(0)) and define the change of energy forward in time 
A NJ (x(l), x(0)) = E N (x f (2),x(l)) - E N (x(l),x(0)), 
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and the change of energy backward in time 

A' NJ {x(l),x(0)) = E N (x(0),x(l)) - E N (x(l),x f (2)). 

Let I = {11,12, ■ ■ ■ ,ik} be the set of all i such that f(xi(l)) ^ g(xi(l)) and fix an 
index j € I. Since e is a conserved quantity for both / and g, we have that 

A' NJ (x(l),x f (2)) - A' Ntg (x(l),x f (2)) = 0, 

and hence we have 



'N-l 



N-l 



J2 €&(()), ^(1)) ~ E <^(1),7(/(^(1)),^(0))) 

i=0 i=Q 
N-l N-l 

e(xi(0),Xi(l)) - 7(5(^(1)), ^(0)) 



i=Q 



i=0 



^ e{xi{l)Mf(*iO))MO)) ~ e(xi(l), 7(5(^(1)), ^i(O))) 
iel\j 

+€(x i (l), 7 (/(x i (l)),x i (O)))-e(x i (l),70/(x i (l)) > x i (O))) 
0. 



(7) 



Consider now x(0) and replace Xj(0) by , y(f(xj(l)),j(g(xj(l)),Xj(0))). We denote 
this new state by x(0) and consider the evolution with initial conditions (x(l),x(0)) 
under the dynamics induced by /. At time t+2, the value of the j th dependent variable 
is equal to 



7(/(%(l)) )7 (/(x i (l)) ) 7(5(%(l)),^(0))))=7(ff(%(l)),^(0)). 



(8) 



Here we have used 7(2/, 7(y, #)) = x, with y = f(xj(l)) and x = ~/(g(xj(l)), Xj (0)). 

By a similar argument, if we consider the evolution with initial conditions (x(l), x(0)) 
under the dynamics induced by g, at time t + 2, the value of the j th dependent variable 
is equal to 

7 (ff(%(l)),7(/(%(l)),7(5(%(l)),^(0)))) 

= 7(^(%(1)),7(5(%(1)),7(/(%(1)),^(0)))) 

= 7(/(x i (l)),x i (0)). (9) 

Here we have used the fact that, by hypothesis, j(x, j(y, z)) = 7(2/, j(x, z)) and there- 
fore 

7 (/(x i (l)),7(ff(%(l))^i(0)))=7(ff(%(l)),7(/(%(l))^i(0))). 

Because of this, for initial conditions (x(l),x(0)), conservation of energy gives rise to 
the identity 

A' NJ (x(l),x(0)) ~ A' N , g (x(l),x(0)) = 0, 



or 



7(/(^(l)), x»(0))) - e(xi(l), 7(5(^(1)), a*(0))) 

+(e(%(l),7(^(l)),^(0)))-e(%(l),7(/(%(l)),^(0))) 
= 0. (10) 
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By subtracting ([7|) from (fTUj) we obtain 

e(%(l), 7 (/(%(l)),x j (0))) = e(%(l),7(3(%(l)),^(0)))- 
Since we can repeat the argument for any j G /, the lemma follows. 



□ 



7 The theorem 

Let us denote by / V g the component- wise maximum of / and g. Then we have: 

Theorem 1 Let f and g be second-order reversible cellular automata of neighborhood 
size (3 and reversible combiner 7, and suppose e = e(xi(t),Xi(t — 1)) is an additive 
conserved density of size a for f and g. Then e is an additive conserved density also 
for h, where h = f V g. 

Proof. Suppose that initial conditions 

x(0) = {x (0),xi(0),...,x^i(0)} 

and 

x(l) = {x (l),xi(l),... ,xat-i(1)} 
are given. Using the notation introduced in the previous section, we need to show that 

A' NJ (x(l),x(0)) = A' NJVg (x(l),x(0)) = 0. 

This amounts to showing that 

N-l N-l 



Y, e{xi{0), Xi (l)) = ^ e (x 4 (l),7((/V 9 )(x,(l)),x. t (0))). 

i=0 i=0 

First notice that we can write the right-hand side of Eq. ljlip as 



(11) 



ieAi 



+ ^e(x i (l), 7 (^(l)),a: i (0))) 



«gA 2 



where 



Ai 
A-2 
A3 



+ X)e(x < (l) >7 (/(£i(l)),x i (0))) J 

igA 3 



{i {mo-)) > g(um 
{i 1/(^(1)) < 5(^(1))} 

{i |/(^(l)) =g(xi(l))}. 



(12) 



By Lemma [3] we have that 

e(S(i,l), 7 (^(S(i,l)),6(i,0))) =e(S(i ) l), 7 (/(S(i ) l)),6(i,0))), 



(13) 
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and therefore the expression in (|12p is equal to: 



£ c(S(i,l),7(/(S(i,l)),6(i,0))) + £ e($(i,l), 7 (/(S(i,l)),&M))) + 




^#,l),7(/(i(i,l)),l(j,0))) = ^e(S(z,l),7(/(^,l)),^0))). 




Since e is a density conserved quantity for / the above is equal to 



N-1 



^e(S(i,0),6(z,l)) 



i=0 



and therefore we have that 



N-1 N-1 



0), 1)) = ^ e(S(i, 1), 7 ((/ V (/)(&(*, 1)), b(i, 0))). 




□ 



Example 2 ^4s i/ie size o/ the alphabet increases, the number of possible reversible 
combiners undergoes a combinatorial explosion, so computing the equivalence set Ti 
becomes computationally impossible. What we can do, however, is to fix a reversible 
combiner and hunt for pairs (/, g) of CA with the same conserved quantities and check 
the validity of the theorem by observing that ip(f) C ip(f V g). 

For instance, for \ S\ =3 and the reversible combiner ^ such that r y(0, 0) = 0; -y(l, 0) 
= 1; 1 (2, 0) = 2; j(0, 1) = 1; 1) = 0; 1 (2, 1) = 2; 1 (0, 2) = 2; 2) = 1; 1 (2, 2) 
= 0, we found that the CA (expressed in Wolfram's notation in base 3) 1311051521973 
and 2889590638889 belong to the same class. Then we computed 1311051521973 V 
2889590638889 = 6589964923167 and verified that ^(1311051521973) C ^(6589964923167). 

Remark 3 The substance of the theorem is that of proving the existence of an order 
relation on CA which is preserved by the conserved quantities. Consider the set J^j of 
second-order reversible CA with a fixed window I with the partial order relation induced 
by V: 



Then, if we consider the set (C) = {/ V g}f j9 ^c generated by any equivalence class 
C G J-j/ ~, the theorem says that if) is an order-preserving map. Indeed, since there is 
a natural partial order C between subspaces ofM. n , we have that 



In this article we have studied a certain class of second-order, reversible CA that were 
first considered by Takesue [1] for their ability to exhibit thermodynamic behavior. 

The inspiration for this work has, in part, come from Noether's theorem, which 
is one of the most celebrated results both in mathematics and physics. It asserts 




f<9^rp(f) ^(j). 



8 Conclusions 
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that conserved quantities are in correspondence with symmetries of the laws of nature. 
Time-translation symmetry is associated with conservation of energy, space-translation 
symmetry with conservation of momentum, and rotation-symmetry is associated with 
conservation of angular momentum. 

These results hold under the assumption that space and time are both continuous: 
to what extent do they still hold for discrete systems such as CA? The answer to this 
question is still an open problem, and a fully satisfactory discrete version of Noether's 
theorem has not yet been formulated. In this work, however, we show that there exists 
a logical organization of additive conserved quantities which reflects some symmetries 
of the system, at least for a large class of second-order, reversible CA. In particular, 
we have shown that there is a partial order relation defined on CA which is preserved 
by the set of their respective conserved quantities. 

There are still many open problems in the theory of reversible CA which we hope 
to approach using our result. In particular, there is an interesting conjecture claimed 
in [2] which states that "antisymmetric conserved quantities", that is those such that 
E(x(t), x(t + 1)) = —E(x(t + l),x(t)), can only be associated to CA endowed with 
special symmetries (Rule 90, for instance, is one of them). More generally, it appears 
that, even in the discrete case, it is possible to classify conserved quantities according 
to symmetry groups. 

As in a previous paper of ours [8], we have herein suggested an approach to this 
problem which only uses combinatorial techniques. We hope that this work will con- 
tribute towards the formulation of a discrete version of Noether's theorem, which so 
neatly classifies the conserved quantities of continuous systems according to the char- 
acteristics of their dynamics. 
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